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We analytically solved the higher-order nonlinear Schrodinger (HNLS) equation with non-Kerr 
nonlinearity under some parametric conditions and investigated explicitly bright and dark solitary 
wave solutions. Periodic wave solutions are also presented. The functional form of the bright and 
dark solitons presented are different from fundamental known sech{.) and tanh{.) respectively. We 
have estimated theoretically the size of the derivative non-Kerr nonlinear coefficients of the HNLS 
equation that agreed the reality of the waveguide made of highly nonlinear optical materials, could 
be used as the model parameters for sub-lOfs pulse propagation. 
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Optical solitons have promising potential to become 
principal information carriers in telecommunication due 
to their capability of propagating long distance without 
attenuation and changing their shapes. Considerable at- 
tentions are being paid theoretically and experimentally 
to analyze the dynamics of optical solitons in optical 
waveguide. The waveguides used in the picosecond op- 
tical pulse propagation in nonlinear optical communica- 
tion systems are usually of Kerr type and consequently 
the dynamics of light pulses are described by nonlinear 
Schrodinger (NLS) family of equations with cubic non- 
linear terms The validity of the NLS equation as a 
reliable model is dependent on the assumption that the 
spacial width of the soliton is much larger than the car- 
rier wavelength. This is equivalent to the condition that 
the width of the soliton frequency spectrum is much less 
than the carrier frequency. The robustness of the optical 
soliton makes it useful for long distance optical communi- 
cation systems, the high frequency of the optical carrier 
makes possible a high bit rate, and to increase the bit rate 
further it is desirable to use shorter femtosecond pulses 
and in order to model the propagation of a femtosecond 
(<100fs) optical pulse in an optical fiber, higher order 
nonlinear Schrodinger equation (HNLS) (not including 
optical fiber loss) [2] 

E, = i{aiEu+a2\EfE)+a3Ettt+a4{\EfE)t+a5Ei\Ef)t 

(1) 

is required. Here z is the normalized distance along 
the fiber, t is the normalized time with the frame of 
the reference moving along the fiber at the group ve- 
locity. The subscripts z and t denotes the spatial and 
temporal partial derivatives respectively. The coefficients 
tti {i = 1, 2, 5), particularly, (ai = 02 = 71, 03 = 
a4 = — ^ and 05 = 71 Tr) are the real parameters 
related to group velocity dispersion (GVD), self phase 
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modulation (SPM), third-order dispersion (TOD), self 
steepening and self-frequency shift due to stimulated Ra- 
man scattering (SRS) respectively. Here (Hj = {^^)u:=ljo 
is the dispersion coefficients evaluated at the carrier fre- 
quency Wo, with the inverse of group velocity, /32, 
the group velocity dispersion parameter, ^3 third order 
dispersion (TOD) parameter and so on. (3 is propaga- 
tion constant. More specifically, 71 is coefficient of cubic 
nonlinearity, which results from the intensity dependent 
refractive index. The term related to — results from the 
intensity dependence of the group velocity and causes 
self steepening and shock formation at the pulse edge. 
The last term related to 05 = 71 Tr incorporates the 
intrapulse Raman scattering and originates from the de- 
layed response, which cause a self-frequency shift, where 
Tr, is called Raman time constant, can be estimated 
from the slope of the Raman gain (SRS). The charac- 
teristic Raman time constant T/j is defined as the first 
moment of the nonlinear response function jSj . Actually, 
to model intrapulse Raman Scattering the last coefficient 
05 should be 05 — iji Tr. However, in the present work 
we succeeded in deriving analytical solutions in the case 
when 05 was real. Let us remark that the latter case 
(a5 real) also dominates in the analytical studies (e.g., 
Painleve property, inverse scattering transform, Hirota 
direct method, conservation laws ) undertaken to date 
Q to show its integrable nature and obtained different 
types of exact solutions such as new solitary wave solu- 
tion, w-shaped solution, bright and dark optical solitary 
wave solutions etc. Thus besides considering that the 
present works Q realizes a significant advance as com- 
pared to previous literature, extension of the work to the 
case of 05 imaginary represents a theoretical challenge 
that should be undertaken in the near future. 

Present days applications in telecommunication and ul- 
trafast signal routiitg systems, as the intensity of the in- 
cident light field becomes stronger, iton-Kerr nonlinear 
effects come into play and due to these additional effects, 
the physical features and the stability of the NLS soliton 
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can change. The way through the non-Kerr nonhnearity 
influences NLS soliton propagation is described by the 
NLS family of equations with higher degree of nonlinear 
terms [sj . To increase the channel handling capacity and 
ultra high speed pulse, it is necessary to transmit soli- 
tary waves at a high bit rate 1 — 10 /s) of ultrashort 
pulses, which can be seen in many applicative contexts 
such as high repetition pulse sources based on fiber tech- 
nology [©]. At the same time, it is also important to in- 
clude some additional higher-order perturbation effects 
to HNLS equation to analyze the solitary wave solution 
in non-Kcrr nonlinear medium. 

Here, in this Letter, we consider the higher-order NLS 
(HNLS) equation with non-Kerr term can be writ- 
ten in terms of slowly varying complex envelope of the 
electric field E{z,t), as 



: i{aiEtt + a2\EfE) + a3Ettt+a4{\EfE)t+a5E{\Ef) 



lae, 



EfE + a7{\EfE)t + asEi\Ef)t 



(2) 



The terms related to coefficients ag, 07, ag in Eq. (2) 
represent the quintic non-Kerr nonlinearities. The quin- 
tic nonlinearities arise from the expansion of the refrac- 
tive index in power of intensity / of the light pulse : 
n = nQ + n2l + n4l'^ + .... Here no is the linear refractive 
index coefficient and n2, are the nonlinear refractive 
index coefficients, originate from third- and fifth-order 
susceptibilities respectively. The polarizations induced 
through these susceptibilities give the cubic and quintic 
(non-Kerr) terms in nonlinear Schrodinger equation re- 
spectively. The nonlinearity arises due to fifth-order sus- 
ceptibility can be obtained in many optical materials such 
as semiconductors, semiconductor doped glasses, Al- 
GaAs, polydiacctylene toluene sulfonate (PTS), chalco- 
genide glasses and some transparent organic materials. 
When the last three terms related to ae, ay, ag of Eq. 
(2) are ignored, the resulting equation becomes the HNLS 
equation as given in Eq. (1). In a recent paper we 
have investigated the Dark-in-the-Bright (DITB) solitary 
wave solution of Eq. (2). The DS or DITB solitary wave 
solution is composed of the product of bright and dark 
solitary waves. We also investigated the stability of the 
DITB solution under some initial perturbation on the 
parametric conditions. The shape of pulse remains un- 
changed up to 20 normalized length even under some very 
small violation in parametric conditions. More recently 
Q , we have studied modulational instability (MI) of Eq. 
(2) with forth-order dispersion in context of optics and 
presented an analytical expression for MI gain to show 
the effects of non-Kerr nonlinearities and higher-order 
dispersions on MI gain spectra. In our study we also 
demonstrate that MI can exist not only for anomalous 
group velocity dispersion (GVD) regime but also in the 
normal GVD regime and also investigated that the quin- 
tic non-Kerr nonlinear terms are more important over the 
cubic Kerr nonlinearity because non-Kerr nonlinearities 
are responsible for stability of localized solutions. But in 
the previous works 

[1,0, [3 the solitary wave solutions in 
presence of higher non-Kerr nonlinearity have not been 



investigated. In a very recent work, Triki and Taha [9| 
presented solitary wave solutions for HNLS equation in- 
cluding non-Kerr nonlinear terms upto the coefficient og 
of Eq. (2). In this article, we shall study the bright- 
and dark solitary wave solutions for HNLS equation that 
contains time derivative of non-Kerr nonlinear terms and 
estimate the size of model coefficients of Eq. (2) which 
will be useful for propagation of very short pulse of width 
around sub-lOfs in highly nonlinear optical fibers. 

To investigate the existence of analytic wave solution 
of HNLS equation in presence of non-Kerr terms we begin 
by scaling the variables of the Eq. (2) in the form 

E = 61* , z = 62^ , and t = b^r 

and choosing 61, 62 and 63 such that the coefficients cor- 
responding to GVD, SPM and TOD become unity. Thus 
in the scaled formed of the Eq. (2) becomes 

*e = j(*rr + 1*1^*) + *rrr + ai(|*|^*)^ + Q2 * ( 1 1 1 ^ )^ 



IQSI*!"** + Q4(|*l'**)r + a5*(|*hr 
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In writing (3) we have chosen bi = (^-^^^ 

To obtain the exact solitary wave solutions of Eq.(3) 
we consider a solution of the following form 

t) = V{t + Oe'^'^?-^^") = 7'(x)e''(^«-^^") , (4) 

with Vix) real. On substitution Eq. (4) into Eq. (3) and 
after removing the exponential terms one can obtain the 
real and imaginary parts of the resulting equation as 



'Pxx = {v-2n + 3fl^)r- 
and 



3Qi + 2Q2 p3 



/2^±+^r^ (4a) 



^ (fc + n^-n^) ^ _ (l^»Qi) ^3 _ (a3-»Q4) ^5 
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In the following we will discuss two cases: 



(4b) 



Case 1: 7^ i. Equating these two equations we 
get the following necessary and sufficient conditions on 
57 and equation for k in terms of fl: 



3ai -|- 2a2 



6{ai 



a2) 



5a4 -I- 4q;5 — 5a3 
(10a4 + 1205) 



(5a) 
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FIG. 1: (Color online) (a)Intensity of the solitary wave profile 
1^(0, r)|^ as a function of r for the value ai — a2 — 1 and 

V = 1 (b) Periodic wave profile |^'(0, t)| as a function of r 
as computed from Eq. (10) for the value ai = 02 = 1 and 

V = -1. 



and 

k = {1 - 3n)(v - 2n - 

with constraint relations 
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The function V{x) satisfies the ordinary nonlinear dif- 
ferential equation in (4a). Multiplying (4a) by 7^^ and 
integrating we get 



(7) 



3ai+2Q;2 
6 



I 5a4+4a5 
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where (v - 2n + 30^ 

and £ is the arbitrary constant of integration. Eq. (7) 
describe the evolution of the anharmonic oscillator with 
potential 



U{V) 



2 2 
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(8) 



and £, the integration constant in (7) is the energy of 
that anharmonic oscillator. For zero energy (£ ~ 0) we 
have find the solution of Eq.(7) as 



\/((a - (fe2 - 4ac) e4v/^x)2 + iQab^ce^^x 



(9) 



provided that a > , 6 > 0. Now using Eqs.(4), (5b) 
and (9) we can write the solution of Eq. (3) as 

^(^ r) = V{t + V^) e<((l-3i2)(t'-20+3J7")-0"+0^)C-Or) 

(10) 

The intensity profile of the solitary wave solution (Eq. 
(10)) is shown in Fig. 1(a), as computed from Eq. (3) 
for the values ai = 02 = 1 and v — 1. One can check the 
evolution of the intensity profile and it is an interesting to 
note that the wave profile remains unchanged during the 
evolution. For negative value of the parameter a, Eq. (3) 
shows the periodic solution. We have presented in Fig. 
1(6) the periodic wave solution. Here we have taken the 
same parameter values as that in Fig. 1(a) but v — —1 
such that a < 0. For = (i.e. 3ai + 2a2 — 3) and 
consequently, a — v, b — ^ and c = ^ one can verify the 
solitary wave solution and periodic solution of Eq. (3) 



FIG. 2: (Color online) Intensity of the solitary wave profile 
|>I'(0, r)| (Eq.(12)) as a function of r as computed from Eq. 
(3) for the value ai — 02 — 1 and = — | 



from Eq. (10) for a > and a < respectively. 
Dark solitary wave. We have seen for £ = 0, HNLS 
equation in presence of non-Kerr terms supports the 
bight optical wave solution provided that the parame- 
ter a in Eq. (7) is positive. Now we chose £ 
c — —j^ such that wc can write Eq.(7) in the form 



ffc and 



dx 



[a-bV^ f 
Sab 



dV 



(11) 



Integrating Eq. (11) we find the solitary wave solution of 
Eq. (3), with the same constraint relations stated above 
in (6), as 



a(T + u^) 



Vb^i^- 



0^ 



(12) 



with fc = (1 - 3n){v ~2n + sn^) + n^, provided 

a < 0, with V = 2n-3n'^ - ^(3q;i + 2a2). In Fig. 3 we 
have plotted the intensity profile of the optical solitary 
wave solution with the parameter values ai = 012 = 1 
and V — We call this solution as a dark solitary 

wave in the sense that the intensity profile associated 
with such soliton exhibits a dip in a uniform background 
and asymptotic absolute value of 4'(^, t) tends toward 
a constant nonzero value for large values of r. In this 
context, although bright solitons are relatively easy 
to generate in optical fiber, the dark solitons are less 
sensitive to optical fiber loss, less infiuenced by noise 
and are more stable against Gordon Haus jitter in long 
communication line [10|. As the mutual interaction 
between two neighboring dark solitons is much weaker 
than that between two bright solitons [ll|, so the 
properties of the dark soliton attracted scientists very 
much in the communication systems. But it is difficult 
to use a dark soliton with a tan/i(. )-type wave form in a 
transmission system because dark pulse cannot be easily 
generated. 

Case 2: f2 = i , Eq. (46) takes the form 
{k + n^ -n^)V -{l-nai)V^ -{a3.-VLai)V^ = . (13) 

Setting the coefficients of 7^, 7''^ and to zero in Eq. 
(13) and using = i we obtain k = —-^ and the con- 
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straint condition 

il = — and ft = — , which implies 04 — aia^ . 

(14) 

In this case, V{x) satisfy the ordinary differential equa- 
tion similar to that in (7) but only differs in the coefficient 
values. 

{V^f ^a'V'^-b'V'^-c'r^ + 2£' . (15) 

Here, in this case, (w - 5) = a', = b', ^oiJ^ ^ 

and £' is the arbitrary constant of integration. In case 
2, similar to Eq. (7) one can solve Eq. (15) to get the 
bright, dark and also periodic wave solutions of Eq. (3). 
Only here the coefficients a', b' and c' are different from 
that in Eq. (7). 

Before conclusion, let us discuss some application of 
the above theoretical prediction. For large channel han- 
dling capacity in the frame of dense time-domain mul- 
tiplexing and for high speed, it is necessary to trans- 
mit solitary wave at a high bit rate of ultrashort pulses. 
At the same time, it is also important to consider the 
higher-order non-Kerr like nonlinearity including deriva- 
tive in HNLS equation for sub-lOfs pulse propagation. 
The relevance of these terms is also important in the 
frame of post-soliton compression that can be achieved 
in highly nonlinear material. Compared to silica glasses, 
as chalcogenide glasses exhibit an extremely high nonlin- 
ear refractive-index coefficient that can be two or three 
orders of magnitude greater than that of silica at 1.55 
/im. They also offer several distinctive optical proper- 
ties such as a transmission window that extends far into 
the infrared (IR) spectral region (up to 25 fim for tel- 
luride glasses). Because of high nonlinearity and large 
IR transparency, chalcogenide fibers are well suited for 
compact Raman amplifiers, supercontinuum generation 
and other mid-IR sources. For experimental verifica- 
tion of the propagation of solitary wave, one may use 
the waveguide made of chalcogenide glasses which are 
made from heavy chalcogen elements such as S, As, Se, 
Te, having electron shells which are easily polarizable 
under an electro-magnetic field excitation. For exam- 
ple, the chalcogenide glass As2Se3, the nonlinear in- 
dex coefficients [H, [l^ are 712 = 2.2 x V^~^'^ rr? jW 
and 714 = —6.5 x \^~^^ jW"^ and for the As^Sz sam- 
ple, the values are = 4.2 x 10~^®to^/M^ and n4 = 
—6.0 X 10~^^ m"^ /W^ . In general, the nonlinear coeffi- 
cients 7i (i = 1, 2) can be estimated from 7^ = ttt^ j 

where rij [j — 2, 4) are the nonlinear refractive index co- 
efficient. Af.ff — TTw'^ is the effective fiber core area, with 
w is the core radius of the fiber, which varies from (3—3.8) 
fim for As2Se3 and (1.3 — 1.7) fim for the AS2S3. A, the 
typical telecommunication wavelength 1.55^m. For the 
range of |7i| = (2000 - 3000)W~^ /Km, cubic nonlin- 
ear coefficient, we have calculated the values for quin- 
tic nonlinear coefficients: I72I — (1.3 — 3)W~''^/Km for 
As2Se3 and for I72I = (3.3- 7.6)W-^/KmAs2S3 re- 
spectively. The 1.55/xm window today is mainly of in- 
terest to long-distance telecommunications application. 



The other important feature of this transmitted wave- 
length is that it matches the fiber's low-loss regions. 
Fiber energy loss (absorption) compensation with suf- 
ficient Raman gain and distortionless propagation of pi- 
cosecond soliton pulses in a monomodal optical fiber have 
been experimentally demonstrated by Mollenauer et. al. 
(l3 |. Chalcogenide glasses have attracted much interest 
in the past few years as a nonlinear optical material in 
the telecommunications wavelength window of 1550 nm, 
and are promising candidates for planar non-linear opti- 
cal (NLO) rib waveguide devices due to high nonlinearity, 
high refractive index, and non-linear optical losses (0.05 
dB/cm) at 1550 n m fisl [T6j . In a vary recent work, M. 
El-Amraoui et. al. [l3| measured the fiber losses as low 
as 0.35 dB/m at 1.55m for a 45 meters long 2.3 m core 
size fiber. The related nonlinear Kerr coefficient is esti- 
mated as high as 2750 W'-^km^-^ . It is also an impor- 
tant to note the fact that chalcogenide glasses exhibit the 
highest nonlinear refractive indices and suffer, at worst, 
only moderately from two-photon absorption, also they 
do not suffer from free-carrier absorption. The nonlin- 
ear absorption faced by the fiber material can be com- 
pensated using derivative higher-order nonlinear Raman 
gain terms. In this context, Tuniz et. al. [l8| studied how 
Raman gain and nonlinear absorption counteract across 
the C and L-bands in two-photon absorption effects in 
single-mode chalcogenide fiber. 

Now, physically, for the ultrashort laser pulse propa- 
gation through optical fiber at telecommunication wave- 
length 1.55 /im, (last reference of Q) and carrier fre- 
quency Wo = 1.22 X lO^^s-i i.e. Tq = 5.1475 x lO^^S, 
if we choose the typical real experimental value for the 
model parameters of Eq. (2) as oi = ^ = 10 ps^/ km, 
a2 = 71 = 2765iy-Vfcm, 03 = = 0.0235 ps^/fcrn, 
a4 = = -14.2328 iy~V((27r)fcTOTi7z) and 05 = 

7iTij = U931W-^ fs/km {Tr = 5.4 /s for chalcogenide 
glass fiber [13|), we can estimate the size of the coef- 
ficients of the non-Kerr nonlinearities of Eq. (2) from 
the constraint relations in Eq.(6). The calculated val- 
ues for the cogf^cients of non-Kerr nonlinearities are 
06 = 72 - ^ = 2.533 iy-VA:m, ay = = 

-1.304 X 10-2 VF-V((27r)fcmriJz) and ag = ^^^^ = 

11.3996 /s/ZcTO. For the sub-lOfs bright pulse com- 
munication in non-Kerr medium with the above model 
parameter values for Eq. (2), one can check that we 
need chalcogenide optical fiber of the type As2Se3 of 
core radius 3.20415/zm and for AS2S3, the core radius 
is 1.399/zm respectively. If we use the same parame- 
ter values for the dark pulse propagation, one can ob- 

2 

tain the energy value for dark pulse will he £ = = 

_ (.-20+3n-)- = -88.779. 

3qi+2q2 

In conclusion, we have reported optical bright and 
dark solitary wave solutions of Higher-Order Nonlinear 
Schrodinger equation in presence of non-Kerr terms 
subject to constraint relations among the parameters. 
The derivative Kerr and non Kerr nonlinear terms are 
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important for compensation of the nonlinear absorption 
during propagation in highly nonlinear material. These 
terms also play an important role for the post-soliton 
compression to get highly stable compressed optical 
pulse. The functional forms of bright and dark solitary 
wave profiles reported are new. We also have presented 
the periodic solutions which are very meaningful in 
optics. We have seen that the estimated values for 
non-Kerr nonlinear model parameters of Eq. (2) agreed 
the reality of the waveguide made of chalcogenide 
glasses. So the inclusion of the non-Kerr nonlinear terms 
in Eq. (1) is justified to describe the sub-lOfs pulse 
propagation in highly nonlinear optical fiber. As well 
as Case 1, we have checked the parametric condition 
in Eq. (14) for Case 2 using the model parameters 
given above. The quintic non-Kerr nonlinear terms in 
contemporary optics become very crucial to the up- 
coming applications in ultrafast signal routing systems, 
double doped optical fiber, optical switching etc. The 



periodic solution can be used to study the formation 
of solitons in the periodic stricture if one consider the 
quintic non-Kerr nonlinearity in fiber Bragg grating 
[19| . and finally, with the calculated parameter values 
for highly nonlinear optical fiber made of chalcogenide 
glasses, the HNLS equation, in presence of non-Kerr 
nonlinear terms as higher-order perturbation, not only 
could find application in broadband telecommunication 
that extends far into the infrared (IR) spectral region for 
the bright and dark optical pulses but also Eq. (2) may 
be a new theoretical model equation for experimental 
designing sub-lOfs optical pulse propagation, which will 
be applicable for the next generation optical fiber using 
chalcogenide type high nonlinear optical glasses. 
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